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Gravitational waves from inspiralling binaries are expected 
to be detected using a data analysis technique known as 
matched filtering. This technique is applicable whenever the 
form of the signal is known accurately. Though we know the 
form of the signal precisely, we will not know a priori its pa- 
rameters. Hence it is essential to filter the raw output through 
a host of search templates each corresponding to different 
values of the parameters. The number of search templates 
needed in detecting the Newtonian waveform characterized 
by three independent parameters is itself several thousands. 
With the inclusion of post-Newtonian corrections the inspiral 
waveform will have four independent parameters and this, it 
was thought, would lead to an increase in the number of fil- 
ters by several orders of magnitude — an unfavorable feature 
since it would drastically slow down data analysis. In this 
paper I show that by a judicious choice of signal parameters 
we can work, even when the first post-Newtonian corrections 
are included, with as many number of parameters as in the 
Newtonian case. In other words I demonstrate that the ef- 
fective dimensionality of the signal parameter space does not 
change when first post-Newtonian corrections are taken into 
account. 

PACS numbers: 04.30. +x, 04.80. +z 



Coalescing systems of compact binaries are the most 
promising sources of gravitational radiation for the 
planned laser interferometric gravitational wave detec- 
tors. As a binary system of stars inspirals, due to radi- 
ation reaction, the gravitational wave sweeps-up in am- 
plitude and frequency. The resulting inspiral waveform 
is often called the chirp waveform. As the wave sweeps- 
up in frequency secular post-Newtonian corrections to 
the phase of the waveform become substantial M . When 
such corrections are incorporated in the templates that 
are used in detecting the gravitational wave signal it may 
be possible to glean useful astrophysical information such 
as the masses of the component stars [[y]|] , their equation 
of state, etc. In addition, observation of several such co- 
alescence events can facilitate an accurate determination 
of cosmological parameters, like the Hubble constant and 
after a few hundred events the density parameter |||}. 
In order to extract information of this kind, it is abso- 
lutely essential that the parameters of the wave form be 
determined very accurately. 

Optimal Weiner filtering or matched filtering is a data 
analysis technique that allows a very precise determina- 
tion of signal parameters g. The method consists in 
correlating the raw output of a detector with a wave- 
form, variously known as a template or a filter, whose 



Fourier transform is the Fourier transform of the signal 
divided by the noise power spectrum. A decision about 
the presence or the absence of a signal is made by look- 
ing at the maximum of the correlation. A template whose 
parameters are exactly matched with those of the signal 
enhances the signal-to-noise ratio (SNR) in proportion to 
the square-root of the number of cycles that the signal 
spends in the detector output, as opposed to the case 
when the shape of the waveform is not known a priori 
and all that can be done is to pre-bandpass filter the 
detector output to the frequency band where the signal 
is assumed to lie, and then look at the SNR for each 
data point in the time domain individually H. For an 
interferometric detector, like the LIGO or the VIRGO, 
operating with a lower frequency cutoff ~ 40 Hz and an 
upper cutoff ~ 1 kHz, this means an amplification in the 
SNR ~ 30-40 for the inspiral waveform from a typical 
binary. This enhancement in the SNR not only increases 
the number of detectable events but, more importantly, 
it also allows a more accurate determination of signal 
parameters — the error in the estimation of a parameter 
being inversely proportional to the SNR. In order to take 
the full advantage of matched filtering it is essential that 
the inspiralling binary waveform, and in particular the 
evolution of its phase, be known to a very high degree of 
accuracy |l[] . A mismatch in the phases of the template 
and the signal can drastically reduce the SNR. Though 
we are confident about the event rate, and hence know 
the probability of the signal being present in the data 
train, we will not know a priori what the signal parame- 
ters are. Consequently, the detector output needs to be 
filtered through a number of templates each correspond- 
ing to a particular set of "test" parameters. The number 
of search templates needed to cover the astrophysically 
relevant range of the parameter space depends primarily 
on the dimensionality of the parameter space. 

Cutler et al. Q pointed out that for the purpose of 
constructing a set of matched filters to detect gravita- 
tional waves from inspiralling binaries it is both neces- 
sary and sufficient to work with the so called restricted 
post-Newtonian gravitational waveform. In this approxi- 
mation one incorporates the post-Newtonian corrections 
only to the phase of the waveform working always with 
just the Newtonian amplitude. Going beyond the re- 
stricted post-Newtonian approximation is not expected 
to change appreciably the magnitude of the statistical 
errors in the parameter extraction, so the restricted post- 
Newtonian approximation can be used to estimate sta- 
tistical errors. However, in the post-detection analysis it 



is necessary to employ more accurate templates since the 
use of just the restricted post-Newtonian waveform would 
give rise to some systematic errors. In the restricted post- 
Newtonian approximation the gravitational waves from a 
binary system of stars, modelled as point masses orbiting 
about each other in a circular orbit, induce a strain h(i) 
at the detector given by 



h(t)=A(irf(t))V 3 cos [<p(t)] 



(1) 



where /(£) - the instantaneous gravitational wave fre- 
quency - is equal to twice the orbital frequency of the bi- 
nary; the constant A involves the distance to the binary, 
its reduced and total mass, and the antenna pattern of 
the detector. The detailed form of A will not be of any 
concern in this paper. The phase of the waveform can be 
schematically written as 



tp(t) = <pN(t) + (fipi N + <ppl-BN + 



(2) 



Here (fN (4) is the dominant Newtonian part of the phase 
and (/3_p«Af represents the nth order post-Newtonian cor- 
rection to it. In the quadrupole approximation we have 
only the Newtonian part of the phase given by M 



ip{t) = tp N (t) 
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Here f(t) is the instantaneous Newtonian gravitational 
wave frequency given implicitly by 



t-t a 
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(4) 



Tpf is a constant having dimensions of time given by 



TN = ^L^-5/3 (7r/ J-8/3 

256 



(5) 



and f a and $ are the instantaneous gravitational wave 
frequency and the phase of the signal, respectively, at 
t = t a . Eventhough it is possible to invert / in terms 
of t we shall continue to work with (||) since it allows a 
straightforward interpretation of the parameter tn and 
later of a similar post-Newtonian parameter. We shall re- 
fer to the time elapsed starting from an epoch when the 
gravitational wave frequency is f a till the epoch when 
it becomes infinite, at which time the two stars would 
theoretically coalesce, as the chirp time of the signal. In 
the quadrupole approximation tn is the chirp time. The 
Newtonian part of the phase, namely equation (0), is es- 
sentially characterized by three parameters: (i) the time 
of arrival t a when the signal first becomes visible in the 
detector, (ii) the phase $ of the signal at the time of ar- 
rival and (iii) the chirp mass M. = (fj^M 2 ) 1 ' 5 , where \i 
and M are the reduced and the total mass of the binary, 
respectively. Note that at this level of approximation the 



phase (as also the amplitude) depends on the masses of 
the two stars only through the above combination of the 
individual masses. Consequently, two coalescing binary 
signals of the same chirp mass but of different sets of 
individual masses would be degenerate and thus exhibit 
exactly the same time evolution. This degeneracy, as we 
shall see below, will be removed when post-Newtonian 
corrections are included in the phase of the waveform. 

How many search templates are needed to cover an 
interesting range of the parameter space if we restrict 
ourselves to the Newtonian waveform? Sathyaprakash 
and Dhurandhar R have made a detailed analysis of this 
question and a typical number of filters they quote is 
about a thousand. They have also pointed out that the 
present-day computer technology is well equipped to fil- 
ter the detector output online. However, as pointed out 
earlier it is not enough to consider just the Newtonian 
waveform. Inclusion of post-Newtonian corrections serve 
dual purpose: On the one hand unless the secular post- 
Newtonian corrections are included in the phase of the 
search templates there would be a severe drop in the 
SNR. On the other hand, and more importantly, in order 
to do interesting astrophysics with gravitational waves it 
is essential to remove the degeneracy in the waveforms 
by including post-Newtonian corrections. 

When post-Newtonian corrections are included the pa- 
rameter space of waveforms acquires an extra dimension. 
It was feared that this would mean a severe burden on 
data analysis: an extra dimension of the parameter space 
implies that it would be necessary to construct for each of 
the thousand odd Newtonian filters a similar number of 
filters corresponding to the post-Newtonian parameter. 
Several authors have therefore analysed the effectiveness 
of a Newtonian template with "wrong" parameters to 
detect a post- Newtonian signal [g]|]] . However, these au- 
thors conclude that even after allowing for a mismatch 
in the parameters of a Newtonian template and a post- 
Newtonian signal the SNR would reduce by about 10-20 
%. Such a drop in the SNR is unfavorable considering 
the low event rate of these sources. In this paper I will 
show that when the first post-Newtonian corrections are 
included in the phase of the waveform it is possible to 
make a judicious choice of the parameters so that the pa- 
rameter space essentially remains only three dimensional. 
While such a strategy is suitable for an unambiguous and 
easy detection it does by no means guarantee a precise 
estimation of all of the binary's parameters. The bot- 
tomline of this paper is that the algorithm presented here 
enables an enhancement in the SNR by including post- 
Newtonian corrections in the search templates without at 
the same time causing any extra burden on data analysis. 
It should however be noted that further off-line analysis 
would be necessary to extract useful astrophysical infor- 
mation. 

With the inclusion of first post- Newtonian correction 
the phase of the waveform becomes p0|,O 



where ipN(t) is given by (g) and 
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Now /(£) is the instantaneous post-Newtonian frequency 
given implicitly by 
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r^r is given by (ph and t p i n is a constant having dimen- 
sions of time given by 



Tpljy 
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(9) 
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The phase un) now contains the reduced mass /i and the 
parameter r\ — fi/M 7 in addition to the chirp mass M.. 
Taking (M., rf) to be the post-Newtonian mass param- 
eters the total mass and the reduced mass are given by 



M = Mrj~ A/b , n = Mri 2/b . Note that the total chirp time 
r of the signal has a Newtonian contribution rjv and a 
post-Newtonian contribution Tpijy : The time left start- 
ing from an epoch when the gravitational wave frequency 
is f a until an epoch when the frequency becomes infinite 
is t = tn + Tpijq. We shall refer to t^ as the Newtonian 
chirp time and to Tpijv as the post- Newtonian chirp time. 
Note that instead of working with the parameters M. and 
r/ we can equivalently work with the parameters tjv and 

Tpl N . 

Thus the post-Newtonian filter is characterized by four 
parameters: Afc = {t a , $, Tjv, Tpijy} where we have used 
the symbol Afe, k — 1, ... ,4, to collectively denote the 
four parameters. Note that t a and $ are kinematical 
parameters that fix the origin of the measurement of time 
and phase, respectively, whereas the Newtonian and the 
post-Newtonian chirp times are dynamical parameters in 
the sense that they decide the evolution of the phase and 
the amplitude of the signal. We shall now set out to see 
if it is possible, for the purpose of filtering, to reduce the 
dimensionality of the parameter space from four to three. 

We begin by defining the scalar product of waveforms 
which plays a crucial role in deciding the filters that are 
required to span the range of parameters and hence to 
assess the effective dimensionality of the parameter space. 
Given two waveforms g(t) and h(t) their scalar product 
is defined by 



(g,h) 



~g(f)h*(f) 

Sn(f) 



df 



(10) 



where S n (f) is the two-sided detector noise power spec- 
tral density and g(f) = f x ' oo g(t)exp(2Trift)dt and 



h(f) = f_ h(i) ex.p(2%ift)dt, are the Fourier transforms 
of the waveforms g(t) and h(t), respectively. The SNR 
p obtained for a signal h(i) using an optimal Weiner fil- 
ter is simply the norm of the signal computed using the 
above definition of the scalar product: 



P= (h,h) 



(11) 



A waveform is said to be normalized if its norm is equal 
to unity. Let us consider the behavior of the scalar 
product of two chirp waveforms g(t;Xk) and h(t;\° k ) 
which differ in all their parameter values, i.e., Afe be- 
ing in general different from A°, and are normalized, i.e., 
{h,h) — (g,g) — 1. Their scalar product C(Afc,A°) is 
given by 



C(\ k ,\° k ) = (g(\ k ),h(\l)) 



(12) 



Here \k can be thought of as the parameters of a signal 
while X° k those of a template. Then C{\ k , \° k ) is the SNR 
obtained using a template that is not necessarily matched 
on to the signal. Since the waveforms are of unit norm 
C(X k , A°) = 1, if A fc = A° and C(X k , A°) < 1, if A fe + A°. 
In general, as indicated by its arguments C(X k ,\1) de- 
pends on the individual values of the parameters both of 
the signal and the template. In what follows I will first 
show that the SNR ( J12| ) depends only on the difference in 
the parameter values Afc — X® . Secondly, I will show that 
a template of a given total chirp time obtains roughly the 
same SNR for all waveforms of the same total chirp time 
though their Newtonian and post-Newtonian chirp times 
may be different from that of the template (see, how- 
ever, the discussion at the end of the paper). The former 
of these two results implies uniformity in the spacing of 
filters [njlj] while the latter result facilitates a massive re- 
duction in the number of templates required in spanning 
the parameter space since instead of constructing filters 
separately for each of the Newtonian and post-Newtonian 
chirp times we can construct filters simply for the total 
chirp time. 

In the stationary phase approximation the Fourier 
transform of the restricted first-post-Newtonian chirp 
waveform for positive frequencies is given by |qj7|,[l3]| 



h(f) = Ar v& 
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For / < the Fourier transform is computed using the 
identity h(—f) — h*(f) obeyed by real functions h(t). 
With the above expression for the Fourier transform the 
SNR Q, using @, takes the form 



C{SX k ) oc 
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where <5Afe = Afc — A?. As in the Newtonian case the 
SNR is independent of the individual parameter values 
of the signal and the template: For all signal-template 
pairs that have the same differences in times of arrival, 
phases, and chirp times one obtains the same SNR. Con- 
sequently, constancy of the distance, measured using the 
scalar product (|1C|), between two nearest neighbour fil- 
ters, that is required in making a choice of filters, trans- 
lates into the constancy of the distance, measured using 
the difference in their parameter values. 

We now seek to analyse the behavior of the SNR Jl4|) 
when the template's parameters are mismatched with 
those of the signal. While it is essential to do the analysis 
for noise power spectral density in real laser interferom- 
eters, such as the one discussed by Finn and Chernoff 
Jl3[ , the results obtained in that case are qualitatively 
the same as in the case of white noise JL4| . In order not to 
divert attention from the main theme of discussion, here I 
will only quote the results for white noise. C(6Xk) traces 
out a four-dimensional surface as we vary 6Xk- In what 
follows I consider the two-dimensional subspace obtained 
by maximizing C(6Xk), over t a and $, for every pair of 
tjv and t p i n of the signal keeping the parameters of the 
template rjy and r plN constant. The surface so obtained 
is plotted for white noise [i.e. S n (f) — const.] in Fig|l|. 
Since the parameters of the template are constant I have 
shown on the x- and y-axis of Fig. [y (and Fig.g) the pa- 
rameters of the signal tjv and t p i n and not the difference 
<5ttv and 6t p i n . The same surface is obtained irrespective 
of what parameters we choose for the parameters of the 
template provided we keep the range of the signal param- 
eters the same. In this sense, the correlation surface in 
Fig. m only depends on the difference in the parameters 
of the signal and the template and not on their abso- 
lute values. For the astrophysically relevant range of the 
masses of the two stars (say, M\,M% G [1, 10] M Q ) and 
for /„ = 100 Hz, t n 6 [4,0.08] s and t p i n e [0.3,0.03] s. 
In Fig.pl the post-Newtonian chirp time is varied over the 
whole of its relevant range while the Newtonian chirp 
time is only varied over a portion of its relevant range. 
The contours of this surface shown in Fig.0 are almost 
straightlines t = tn + t p in = const., except for the 
inner most one or two contours. The value of C corre- 
sponding to the inner most contour is 0.9 and reduces by 
0.1 with successive outer contours. 

A useful interpretation of the surface in Fig.fi] is the 
following: Imagine that we have a template with parame- 
ters corresponding to the center of the grid. It obtains an 



SNR of unity with a signal whose parameters are exactly 
matched on to it. The SNR that it obtains for other wave- 
forms is in general less than unity but, as is evident from 
Fig.p], the SNR is almost equal to unity for every wave- 
form whose total chirp time is the same as its own total 
chirp time. Thus, if we choose our search templates along 
the curve perpendicular to the contours of Fig.0 then we 
will in effect be covering the entire subspace of the sig- 
nal. In other words we can span the two-dimensional 
(rAT,Tpijv) subspace of the four- dimensional parameter 
space with just one parameter. This curve is an appropri- 
ate one of a family of straight lines tn — t p i n + const. 
Consequently, as far as the choice of filters is consid- 
ered we need only work with three parameters, namely 
{t a ,&,r}. This reduction in the effective number of pa- 
rameters can be traced to the fact that there is a strong 
covariance between the parameters tn and t p i^. This 
result, together with the details of a Monte Carlo simu- 
lation demonstrating the effectiveness of the claim made 
in this paper and the results of including higher-order 
post-Newtonian corrections, and noise envisaged in real 
interferometers will be published elsewhere |L4fl . 

Let me conclude by making two cautionary remarks. 
The first one concerns the choice of parameters: It should 
not be thought that the effective dimensionality of the pa- 
rameter space is three only when the set {t a , $, tn, t p in} 
is employed in constructing a lattice of filters. Afterall 
the reduction of dimensionality is related to the property 
of the scalar product ( |l0| ) which is re-parametrization in- 
variant. The advantage of the set used in this paper is 
that it allows us to conclude about the effective dimen- 
sionality without recourse to complicated mathematical 
analysis. However, the final justification has to come 
from a more rigorous analysis which will be taken up 
in a future paper 14 1. The second comment is about 
the scope of the reduced dimensionality of the param- 
eter space itself: The parameter space would be truely 
three-dimensional provided the correlation function is a 
constant in the direction tjv + t p in = const. However, 
as can be inferred from the Fig. [I] and g, strictly speak- 
ing, this is not the case: The correlation function slowly 
decreases as we move away from the maximum of the 
correlation function in the direction tjv + t p i n = const. 
This means two things: (i) The argument about the re- 
duction in dimensionality is only valid as long as the cor- 
relation function has not dropped too much (ideally less 
than about 1%) along lines of constant tn + t p in and 
(ii) A post-Newtonian filter of a given total chirp time 
cannot be replaced by a Newtonian filter of the same 
chirp time. In other words, the presence of the post- 
Newtonian term cannot be mimicked by a Newtonian fil- 
ter alone. However, for the astrophysically relevant range 
of the parameter t p i^ it turns out that we need only use 
a three-dimensional lattice of filters or at worst two sets 
of a three-dimensional lattice. 

I am indebted to the referee for pointing out several 



flaws in an earlier version of this paper. It is a pleasure FIG ' 2 ' Contours of the SNR surface shown in Fig.l. The 

to thank the members of the gravitational wave group at contou rs, are approximately straightlines r N + r plN = const. 
IUCAA, especially R. Balasubramanian, Sanjeev Dhu- 
randhar and Kanti Jotania, for many useful discussions. 
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FIG. 1. Surface showing the maximum, over t a and $, of 
the SNR C(S\ k ). 
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